Unconventional superconducting gap in underdoped cuprates 
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A generic theory of the quasi-particle superconducting gap in underdoped cuprates is derived in 
the strong coupling limit, and found to describe extremely well the experimental "second gap" in 
absolute scale. In drastic contrast to the standard theories of Bogoliubov quasi-particle excitations, 
the quasi-particle gap is shown to originate from anomalous kinetic process, completely unrelated 
to the pairing strength. Furthermore, the fc-dependence of the gap deviates significantly from 
the pure d^2_y2-wave of the order parameter. Our study reveals a new paradigm for the nature of 
superconducting gap, and is expected to reconcile numerous apparent contradictions among existing 
experiments toward a more coherent understanding of high-temperature superconductivity. 
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Recent exciting discovery of the superconducting gap 
(SCG) in the underdoped cuprates by angle-resolved 
potoemission spectroscopy (ARPES) [THU] and scan- 
ning tunneling microscopy (STM) P^BHTO] reveals critical 
clues for the puzzling high-temperature superconductiv- 
ity, which has proven to be one of the most important 
yet challenging problems of condensed matter physics for 
more than two decades. Indeed, unlike the larger pseudo- 
gap (PG), the SCG closes exactly at the transition tem- 
perature Tc [12] , and shows strong correlation to the dop- 
ing dependence of Tc |31I11[2DJ[5T]. Yet, some experiments 
found intriguing indications that the fc-dependence of the 
SGG might not follow the well-established pure d-wave 
of the order parameter |12| . Obviously, a deeper un- 
derstanding of the properties of the SCG in the quasi- 
particle (QP) excitation spectra holds the essential key 
to a resolution to the long-standing problem of high- 
temperature superconductivity in underdoped cuprates. 

Surprisingly, despite the intensive experimental stud- 
ies, to date there have been limited attempts [22l - [26] to 
understand the SCG, in great contrast to the numerous 
efforts to address the PG phenomenon (571 HH] ■ Further- 
more, the existing understanding for the SCG remains in 
the scope of the weak coupling BCS framework, which 
was proposed insufhcient [29] in the underdoped regime 
to account for the essential phase fluctuation of the order 
parameter. As a result, most of the essential questions 
remain open, including the energy scale that controls the 
size of SCG, the true nature of the sharp QP at the edge 
of the SCG, and the precise fc-dependence of the SCG, 
among others. 

In this letter, we address these key questions of the 
SCG by deriving rigorously a generic description in the 
strong coupling limit, in which real-space pairs of holes 
are assumed tightly bound in nearest neighboring sites. 
With a rigorous separation of Hilbert space into bosonic 
and fermionic portions that describe the bound pairs 
and unbound doped holes, respectively, the low-energy 
fermionic coupling to the Bose condensate is made ex- 



plicit. The resulting size of the SCG is shown to describe 
accurately the experimental gap in absolute scale with- 
out any free parameter. Intriguingly, completely different 
from the standard Bogoliubov QP excitation, the SCG is 
found to originate from the anomalous kinetic process 
unrelated to the pairing strength. Furthermore, the k- 
dependence of the SCG deviates significantly from the 
pure (i^2_j,2-wave of the underlying order parameter. The 
new paradigm fills a long-standing void in our knowledge 
on the SCG in the opposite limit of the BCS theory, and 
is expected to reconcile existing and future experiments 
on the QP SCG toward a comprehensive understanding 
of high- Tc superconductivity. 

Since the SCG occurs at lower energy and lower tem- 
perature than the PG physics, we proceed as following. 
First, the average kinetics of the fully dressed one-particle 
propagator are obtained from fitting the experimental 
measurement of the spectral functions in the "normal 
state" PG phase where coherent pairs are not yet formed. 
Second, the many-body Hilbert space of the doped holes 
is split rigorously into those spanned by paired holes 
(bosons) and the unpaired holes (fermions). Finally, the 
QP SCG is obtained from the fermion's coupling to the 
coherent pairs formed below Tc. 

Table I shows the tight-binding parameter obtained 
from fitting the available ARPES data [30l[31] in the PG 
phase of La2_A-Sr5Cu04 over a large energy ('^O.TeV) 
range. The resulting Hamiltonian, H = '^nitnic^^Cii , 
should be properly understood as a convenient represen- 
tation of the averaged one-particle kinetic process given 



5{%) 5.2 7 15 22 

t -4.0 -3.7 1.1 8.6 

t' 29.8 30.6 33.5 35.1 

t" 29.8 29.0 25.8 23.9 

TABLE I: Doping dependent hopping parameters (in meV) 
extracted from ARPES in the PG phase of LSCO [501 l3T] . 



by the dressed one-particle propagator under the full 
renormalization of the many-body interactions. It cer- 
tainly does not provide the explicit information of the 
pairing interaction, nor the decay of the QP. As clearly 
shown in Table I, the fully dressed band structure is dop- 
ing dependent, reflecting the competition between the 
kinetic and potential energy. The nearest neighbor hop- 
ping t is almost entirely suppressed, understandable from 
the strong anti-ferromagnetic correlation, and is thus 
dropped from now on. Interestingly, the value of the 
second and the third neighbor hopping, t' and t" , ap- 
proaches the same near the 5.2% quantum critical point 
(QCP) where superconductivity ceases to exist even at 
temperature T = 0. The equivalence of t' and t" turns 
out to have profound influences on the behavior of the 
superfluid (c.f: Ref [5^) and the SCG as shown below. 

Next, consider the strong coupling regime where the 
binding of pairs are stronger than the fully renormal- 
ized kinetic energy, such that at low temperature and 
low energy, the paired holes cannot break apart easily. 
Such regime can result from the above suppression of 
the kinetic process or the large potential energy associ- 
ated with the anti-ferromagnetic correlation [33] and/or 
the formation of bi-polarons [53] . (As shown below, the 
detail of pairing mechanism is actually not essential for 
the low-energy physics of SCG in this case.) Consistent 
with the pure d^2„j^2 symmetry of the order parame- 
ter, in such a strong coupling regime, the paired holes 
are expected to remain primarily first-neighbor to each 
other as real-space pairs, since doubly occupying a site 
by doped holes is unlikely in a lightly doped Mott insu- 
lator. Similarly, the probability of finding unpaired holes 
below or around T^ is small. It is thus convenient to split 
the many-body Hilbert space into two components: one 
spanned by the tightly bound pairs of nearest neighbor- 
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ing holes, &|- = cLc:-, , and the other by the unpaired 
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holes, fl^. For simplicity, a strong anti- ferromagnetic 
correlation between the nearest neighboring doped holes 
is assumed here (well justified in the underdoped regime). 

Obviously, the two-site nature of the bosons introduces 
potential ambiguities in the representation. To ensure 
a unique representation, the following two conventions 
are introduced: (A) whenever possible, bosonic repre- 
sentation is employed prior to the fermionic one (c.f.: 
Fig. [ija)), and (B) if multiple representations still exist, 
the one that maximizes the total "preference factor" is 
employed, where the preference factor is assigned to each 
pair according to the convention defined in Fig. WCb) . 

In the new representation, the original creation opera- 
tor of a doped hole now corresponds rigorously to 
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FIG. 1: (Color online) Illustration of (a) convention (A) and 
(b) preference factor for a unique representation. (c)(d) Illus- 
tration of two possible contributions of hole creations. 
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where ^ denotes "sum only the first non-zero" accord- 
ing to the convention defined above. Figure lljc) (d) il- 
lustrates the meaning of this non-standard expression. 
When creating an additional hole, one can either (first 
term) create an unpaired hole if and only if no other 
unpaired hole is next to it (Fig. flTc)), or (second term) 
meet an unpaired hole and form a new pair (Fig. fUd)). 
It is important to emphasize the mutual exclusion of the 
terms in Eq. (II]) according to the above conventions [35] . 
Owing to the low-density of the doped holes, the higher 
order terms can be dropped safely. For a cleaner presen- 
tation, the spin indices and the associated sign changes 
will be suppressed from now on, and addition of hermi- 
tian conjugate will be implied. 

The effective one-particle kinetics then turns into: 

H = H^ + Hf + H^f 

= C£ij,A,b^'i) + '^U,f}h') (2) 

ii' j ii' 

ii' j ii' jj' 

where H'' describes the pivoting motion of the pairs 
that results in the local d-wave symmetry [32], H^ 
gives the motion of unpaired holes, and H^^ describes 
the coupling between the bosons and the fermions. 
Obviously, at T > Tc the one-body propagator 
G(i,0) ^ (c(t)ct(O)) = (/(t)/t(0)) -t- (/t(t)6(i)/t(0)) + 
{fit)bHO)f{0)) + (/t(i)6(i)6t(0)/(0)) recovers the orig- 
inal experimental dispersion by construction. (Here () 
denotes time-ordered ensemble average.) At low tem- 
perature, the last three terms of G are insignificant at 
low-energy considering 1) the large pairing/depairing en- 
ergy they need, and 2) the insufficient supply of unpaired 
holes. This leaves only the first term in G essential to the 
low-energy physics of the QP, corresponding to injection 
of an unpaired hole followed by removal of the same at 
some other time. 



At T < Tc, the formation of superfluid of coher- 
ent pairs introduces additional effects to the low-energy 
fermionic Hilbert space through the t^^ terms. Indeed, 
converting this term to momentum space with a dx2_y2 

phase of the order parameter {b\ j,/^_j,), and take the av- 
erage of X]' -> i Z]> one obtains V^Sfe ^V fkf-k + c.c. 
(where hq denotes the number of pairs per Cu atom in 
the condensate and can be practially approximated by 
the superfluid density Ug in mimber/Cu) with 
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{cos{kx) — cos{ky)) 
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X [4t'cos(fca;) cos{ky) + 2t" {cos{2kx) + cos{2ky))]. 

giving the fc-dependence of the coupling. The standard 
Bogoliubov transformation then leads to our main result 
on the low-energy QP gap: 



A.(r) = 4^ 



VMT). 
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Several striking features emerge from this rigorous and 
generic derivation of the QP SCG in the strong cou- 
pling regime. First, in drastic contrast to the weak 
coupling BCS theory, the SCG is entirely controlled by 
the strength of the anomalous kinetic process, t''-^ , and 
completely unrelated to the pairing strength! To ver- 
ify this surprising conclusion, we compare in Fig. K[a) 
our result directly with the ARPES [TT] experimental 
SCG of La2_5Sr5Cu04 (LSCO) at T = 7.6K - 0.54Tc 
(Tc = UK) for 5-7%. The theoretical n^iT) is es- 
timated from ^71, where the thermal reduction factor 
d{T) = n,{T - 0.54Tc)/ns(0) - 0.72 is directly taken 
from the penetration depth measurement of YBCO at 
the same doping level [57] . (We are unaware of similar 
measurement for LSCO.) The quantum reduction factor 
7 accounts for the drastic suppression of the superfluid 
density in the vicinity of the 5.2% QCP, but away from 
the QCP it should approaches 1 from below rapidly. Here 
7 ~ 0.85 is taken from Ref [32] for <5 = 7%. 

Clearly, the direct comparison in Figure [2] shows an ex- 
cellent agreement with the experimental SCG in absolute 
scale. Considering that not even a single free parameter 
is needed in the theory, this degree of agreement in the 
size of SCG is truly remarkable. Evidently, the SCG 
measured by ARPES and STM (and indirectly by other 
inelastic measurements) indeed reflects the strength of 
the effective kinetics of the doped holes, but is completely 
unrelated to the pairing strength, a feature unimaginable 
from the weak coupling BCS theory. 

The 5- and fc-dependence of the SCG holds another 
surprise in our result. As already hinted near the end 
of the experimental data set {4> — 30° ) in Fig. l2ja) , the 
fc-dependence of the theoretical SCG deviates dramati- 
cally from that of the underlying pure dx2_y2 symmetry 
of the order parameter (presented as a fit via the green 
line) . Unexpectedly, the size of the SCG actually reduces 
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FIG. 2: (Color online) (a) Momentum dependence of the SCG 
for 7% doped LSCO at 7.6K. (b) Theoretical fc-dependence 
of the SCG from overdoped (OD/22%), optimally doped 
(OP/15%) to underdoped (UD/5.2%) systems, (c) ARPES 
normalized weight transfer [AW oc A^) in Bi2201. Blue and 
pink arrows indicate the trends upon reduced doping. 



near the anti-nodal direction, instead of reaching its max- 
imum. A more clear picture is demonstrated in Fig. ^b) 
for three doping levels, 5 — 22% (overdoped), 15% (opti- 
mally doped), and 5.2% (extremely underdoped). For a 
clearer illustration of the trend, ng is simply set to 5/2. 
(Thus, an exaggeration of the overall scale of the SCG 
for the overdoped case is to be understood.) 

Two apparent trends can be easily identified from 
Fig. ^h). First, the relative size of the SCG at (7r,0) 
decreases with reduced doping level. In fact, Eq.[3| indi- 
cates that near the QCP at J — 5.2%, where t' and t" 
approaches the same value (c.f.: Table I), the size of the 
SCG reduces to exactly zero at (7r,0), generating addi- 
tional nodal points. Second, the fc-point where the SCG 
reaches its maximum on the Fermi surface approaches 
(7r/2,7r/2) as the doping level decreases, in a manner sim- 
ilar to the end of the Fermi arc. 

Amazingly, both trends are observed in the re- 
cent ARPES analysis [T^ of the normalized spectral 
weight transfer AM/ from T > T^ to T < T^ in 
(Bi,Pb)2(Sr,La)2Cu06+5, shown in Fig.[2];c). (The prop- 
erly normalized spectral weight transfer is roughly pro- 
portional to the size of the SCG and is likely a much more 
reliable measure of the SCG than an artificial assignment 
of a gap edge, given the limited experimental resolution 
and the very low residual spectral weight in the large PC 
near (7r,0).) Indeed, other than the expected exaggera- 
tion of the overall magnitude of the over-doped case, our 
theoretical results resemble very much the experimental 
((5,fc)-dependence, especially capturing the two apparent 
trends unexpected from a pure d^2_y2 symmetry of the 
underlying order parameter, another feature inconceiv- 
able from the weak coupling BCS theory. 

These "anti-BCS" features found in our resulting SCG 
only reflect a new paradigm of the QP SCG in the strong 
coupling limit, as revealed in our derivation. Indeed, 
the smaller SCG observed inside the PC has a novel 
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FIG. 3: (Color online) (a) An example of anomalous cou- 
pling via pairing interaction V in the BCS theory, with green 
lines denoting G. (b) Illustration of Bogoliubov QP. (c) The 
anomalous kinetic coupling in this work, with thick blue line 
denoting bosonic propagator, and thin blue ring emphasizing 
the anti-hole's strong binding to the bosonic space, (d) Illus- 
tration of low-energy QP not involving pairing energy scale. 



cuprates. The resulting SCG is shown to differ funda- 
mentally from the standard Bogoliubov QP excitation in 
the BCS theory. Not only its size is controlled by the 
anomalous kinetic process completely unrelated to the 
pairing strength, but also its ((S,/c)-dependence deviates 
significantly from the d2^2_j^2-wave of the underlying or- 
der parameter. The new paradigm fills the long-standing 
void in our knowledge on the SCG in the opposite limit 
of the BCS theory, and reveals essential nature of the 
high-Tc superconductivity. 

We thank A. Kaminsky for providing the experimental 
data in Fig. Wic). This work was supported by the U.S. 
Department of Energy, Ofhce of Basic Energy Science, 
under Contract No. DE-AC02-98CH10886. 



character entirely different from the well-understood Bo- 
goliubov QP excitation in weak coupling theories (c.f.: 
Fig. [3]ja)(b)) that breaks apart a particle from the pair 
with an energy scale of the binding energy. The coupling 
Aj^^fjfi' = 4/(5i^y )Vj' allows also another low-energy 
anomalous kinetic process (Fig. [3]jc)(d)) of an unpaired 
hole /^ hopping and morphing into an anti-hole g^ = f 
inside the Bose condensate (BEC), and its hermitian con- 
jugate allows the reverse process of anti-hole hopping into 
a hole outside the BEC. Since the creation of the com- 
posite object {gjbijY does not involves the energy of the 
pairing strength (unlike the creation of a single b^), the 
whole effective scattering process (Fig. f3]jc)(d)) only re- 
flects the energy scale of the kinetics t''^ 

An important implication of the above new paradigm 
is that the QPs around the SCG are distinct from the 
holes that form the pairs. Thus, they cannot participate 
the decay of pairs via the standard depairing process [5B] 
advoated previously to explain the near linear reduction 
of the superfluid density against temperature [37] . This 
leaves thermal phase fluctuation |29] or thermal deple- 
tion of the BEC 'S? the only known explanation. Fur- 
thermore, even though both Tc [32 and SCG are un- 
related to the pairing strength, they are correlated via 
n,(0): Tc oc n,(0) [29l [32] and A ex ^/n,{0), again defy- 
ing A/Tc - 1.8 from the BCS theory. 

Obviously, the novel nature of the SCG suggests 
strongly revisiting the existing theoretical considerations 
and experimental interpretations of the SCG. In particu- 
lar, the decoupling of the fc-dependence of SCG from that 
of the order parameter raises an alarming flag on the com- 
mon practice of equating the nodes in the SCG to those 
of the order parameter. Also, the above generic consid- 
erations can be applied to the analysis of other strong 
interaction theories j351 [SH] to provide further physical 
insights unapparent from their numerical solutions. 

In conclusion, we develop a generic theory of SCG 
in the strong coupling limit suitable for underdoped 
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